AD-A035  703 


unclassified 


VANOERBILT  UNIV  NASHVILLE  TENN  DEPT  OF  CHEMISTRY  F/G  20/10 

APPLIED  QUANTUM  CHEMISTRY  OF  NONMETALLIC  MATERIALS. (U)  ' 

DEC  76  J R VAN  WAZER  AF-AF0SR-226S-72 

AFOSR-TR-77-0063  NL 


r 

■rfr; 

PH 

n 

n 

n 

n 

n 

n 

n 

n 

UJ 

U 

LJ 

LJ 



LJ 

LJ 

LJ 

LJ 

LJ 

ADA035703 


FINAL  SCIENTIFIC  REP&T. 


Covering  the  Period  from  February  1,  1972  Through  October  31,  1976 


APPLIED  QUANTUM  CHEMISTRY  OF  NONMETALLIC  MATERIALS 


Ap5.1c.va;]  fov  public  roleeae  j 
cli  strib ut i on  uni lai t ©d « 


Submitted  to  the 


UNITED  STATES  AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 

5-73^ 


VandeAbttt.  UniveAAtty  — 
Ncuhvttle,  Tennessee  37235 


pi  miipanum  u^ator 

Professor  of  Chemistry 
Box  1521,  Station  B 
Vanderbilt  University 
Nashville,  Tennessee  37235 
Telephone:  (615)  322-2790 


Ss> 


sclIS,mc  tusw 

S&5S  S iu“uoa‘ 


r\ £\ 


DepartMnt 

Vanderbilt 


AF  Office  of  8c 
Bolling  A*B,  B1 


Washington,  DC 


it.  supplcmcmtary  JEtsI 


•>.  nr  « j 


»>• 


/i 


9338-02,  61102*^681303 


Approved  for  public  releeee;  distribution  unlimited 


Start  M.  (f  Stftarant  hem  Srtof O 


17.  OllTmaUTION  STATSMSMT  taf  *• 


riTTnn:; 


L J 


f 


r 


I. 

II. 

III. 

IV. 


TABLE  OF  CONTENTS 


ABSTRACT 1 

INTRODUCTION  2 

BRIEF  REVIEW  OF  THE  WORK  PUBLISHED  UNDER  THIS  GRANT 3 

A.  Ab  Initio  Studies  on  Various  Molecules,  With  All  Core 

Electrons  Included  3 

B.  Calculations  with  the  Core  Electrons  Omitted  6 

REFERENCES  CITED  13 

APPENDIX 15 


ABSTRACT 


The  major  advance  made  in  the  work  supported  under  this  Grant  has 
been  the  development  of  what  appears  to  be  a very  generally  applicable 
procedure  for  carrying  out  ab  initio  calculations  employing  only  the 
valence  orbitals  of  atoms,  molecules  and  presumably  also  of  solids, 
and  using  a potential  function  to  describe  the  charge  distributions 
of  the  cores  and  the  Phillips-Kleinman  pseudopotential  technique  to 
account  for  the  other  interactions  between  the  fully  described 
valence  orbitals  and  the  emulated  core  orbitals.  This  new  method  is 
now  being  subjected  to  further  extension  and  development  under  our 
new  continuation  Grant  AFOSR-77-3145.  During  the  period  of  Grant  AFOSR- 
72-2265,  some  advances  were  also  made  in  describing  quantum-chemical 
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APPLIED  QUANTUM  CHEMISTRY  OF  NONMETALLIC  MATERIALS 
I . INTRODUCTION 

The  underlying  .physical  lews  necessary  for  the  mathematical 
theory  of  a large  part  of  physics  and  the  whole  of  chemistry  are 
thus  completely  known  and  the  difficulty  is  only  that  the  exact 
application  of  these  laws  leads  to  equations  much  too  complicated 
to  be  soluble.  It  therefore  becomes  desirable  that  approximate 
practical  methods  of  applying  quantum  mechanics  should  be 
developed  which  can  lead  to  an  explanation  of  the  main  features 
of  complex  atomic  systems  without  too  much  computation . 

- P.  A.  M.  V-Oulc,  Paoc.  Roy.  Soc.,  123 , 714  (1929). 

In  the  Grant  Proposal  for  AFOSR-72-2265,  the  proposed  studies 
were  divided  into  the  following  groups:  (1)  improved  SCF  calculations, 

(2)  computations  on  various  molecules,  (3)  development  of  a priori  methods, 
and  (4)  calculation  of  properties.  During  the  period  of  this  Grant,  some 
effort  was  expended  on  all  of  these  areas,  with  emphasis  on  items  1 and  2. 
At  the  time  this  Proposal  was  written,  it  was  very  clear  to  us  that  a 
formidable  major  problem  in  applying  reasonably  reliable  quantum  mechanical 
calculations  to  all  of  the  various  elements  of  the  Periodic  Table  lay  in 
the  fact  that  such  a large  proportion  of  the  total  computational  effort 
must  be  devoted  to  optimizing  the  core  electrons,  even  though  the 
mathematical  description  for  the  core  of  a given  atom  varies  only  slightly 
from  one  molecular  environment  to  another.  The  reason  for  this  is  that 
quantum  calculations  are  usually  carried  out  by  energy  optimization  and 
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the  core  electrons  make  by  far  the  largest  contribution  to  the  total 


electronic  energy.  For  example,  even  in  the  case  of  a Second-Period 
element  (such  as  carbon)  which  has  only  a single  pair  of  core  electrons 
per  atom,  about  75%  of  the  cost  of  a full  self-consistent-field  (SCF) 
molecular  calculation  is  due  to  the  presence  of  the  core.  For  the 
heavy  atoms,  the  high  cost  of  a full -SCF  study  is  nearly  all  due  to  the 
presence  of  the  core  electrons. 

In  the  original  proposal,  we  suggested  several  ways  for  avoiding 
the  core  problem.  These  suggestions  lay  in  two  classes:  one  consisting 

of  some  way  of  modeling  the  core  electrons;  and  the  other,  some  improvement 
in  the  semi-empirical  approach  whereby  the  core  is  simply  ignored  by  ad- 
justing the  nuclear  charge  so  as  to  correspond  to  the  sum  of  the  charge 
of  the  nucleus  and  of  the  omitted  core  electrons. 

II.  A BRIEF  REVIEW  OF  THE  WORK  PUBLISHED  UNDER  THIS  GRANT 

A.  Ab-Initio  Studies  on  Various  Molecules t With  all  Core  Electrons 
Included.  The  first  six  papers  supported  under  this  Grant  dealt  with 
ab  initio  calculations,  using  a moderately  sized  Gaussian  basis  set,  on 
a collection  of  molecules  chosen  because  of  their  special  chemical  interest. 

In  our  study  of  the  PFaH  and  OPF,H  molecules,1  it  was  found  that  9 of  the  10 
valence-shell  molecular  orbitals  of  PFaH  are  individually  closely  related  to  9 
of  the  12  valence-shell  molecular  orbitals  of  0PFaH,  with  this  close 
relationship  showing  up  in  the  electronic  population,  electron— density  plots, 
and  orbital  energies.  In  addition,  the  valence-ehell  molecular  orbitals 
were  correlated  with  those  of  the  PF,,  OPF,,  PH,,  and  OPH,  molecules;  and 
it  was  shown  that  some  of  the  molecular  orbitals  persist  practically  unchanged 
from  one  molecule  to  another  in  this  set  of  not  closely  related  molecules. 

In  the  study  of  the  electronic  structures  of  phosphirane  and  thiirane,* 
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1 

it  was  shown  that  the  valence  orbitals  of  the  cyclic  molecules,  C2H<,PH  and 
CaH*S  are  closely  related  and  that  these  are  quite  similar  to  the  respective 
orbitals  of  cyclopropane.  These  two  papers 1,3  indicate  quite  clearly  that 
not  only  are  the  delocalized  molecular  orbitals  obtained  from  self-consistent- 
field  calculations  about  as  readily  understandable  and  interpretable  in 
chemical  terms  as  are  the  localized  orbitals  derived  from  them,  but  that 
some  of  the  delocalized  orbitals  of  one  molecule  may  reappear  practically 
unchanged  in  a related  structure.  We  concluded  in  these  studies  that 

1 

there  are  probably  a relatively  limited  number  of  SCF  delocalized  orbitals 
which  are  to  be  found  over  and  over  again  in  common  chemical  compounds. 

It  was  suggested  that  these  "molecularly  invariant"  orbitals  should  be 
of  great  importance  in  practical  chemistry  and  that  their  contributions 
to  properties  (such  as  the  electron-binding  energy,  the  contribution  to 
the  electronic  part  of  the  dipole  moment,  etc.)  ought  to  be  intercompared 

1 

to  give  a better  insight  into  the  electronic  structure  of  matter. 

I 

In  an  ah  inito  study3  of  the  role  of  d orbitals  in  chlorosilane,  the 
electron-density  plots  demonstrated  that  conferring  d character  upon  either 
silicon  or  chlorine  has  about  the  same  effect  on  the  detailed  transfer  of 
charge  between  the  chlorine  and  silicon  atoms.  However,  the  usual  population 
analysis  leads  to  the  conclusion  that,  when  d character  is  allowed  only 
to  the  silicon,  there  is  a p -d  transfer  of  charge  from  the  chlorine  to 

n n 

the  silicon  but,  when  it  is  allowed  only  to  the  chlorine,  the  effect  is 
essentially  Just  a polarization  of  this  atom.  In  other  words  for  this 
example  (and  presumably  in  the  usual  case  ),  there  is  really  essentially  no 
physical  difference  between  Pff~’dr  charge  transfer  from  atom  B to  atom  A 
and  simple  polarization  of  the  charge  of  atom  B towards  A. 

It  was  also  pointed  out  in  this  work  that  the  calculated  variations 
in  the  orbital  enarglas  of  the  core  electrons  upon  allowing  or  disallowing 
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d character  to  the  silicon  or  chlorine  atoms  may  be  interpreted  in  terms 
of  changes  in  the  electrostatic  potential  in  the  core  region.  These  findings 
have  significance  with  respect  to  the  Interpretation  of  data  from  inner- 
orbital  photoelectron  spectroscopy. 

The  next  two  papers4'9  deal  with  the  role  of  d orbitals  in  a sulfur 
and  a phosphorus  compound,  employing  ab  initio  SCF  calculations  with  a moderately 
sized  Gaussian  basis  set.  In  the  first  of  these  studies,4  the  electronic 
structure  of  isothiocyanic  acid  was  calculated  and  then  compared  with  that 
calculated  for  lsocyanic  acid;  and  it  was  concluded  that  (a)  the  it  system 
in  HNCS  involves  a nitrogen  lone  pair  stabilized  by  a higher-lying  C-S  it 
bond,  while  the  r.  system  of  HNCO  consists  of  a C-0  tt  bond  stabilized  by 
the  higher-energy  nitrogen  lone  pair  and  (b)  the  d orbitals  of  sulfur  are 
accepting  density  in  a o rather  than  a ir  fashion.  Although  we  put  a tremendous 
amount  of  effort  in  our  SCF  study  of  the  electronic  structures9  of  phosphorus 
pentafluoride  and  tetrafluorophosphorane,  it  turned  out  that  we  were  scooped  in 
our  publication  by  Strich  and  Veillard.*  Since  these  authors  employed  a larger 
basis  set  than  we  did,  our  publication9  as  it  finally  appeared  represented 
a considerable  condensation  of  our  full  study.  Not  surprisingly,  our 
calculated  results  were  quite  similar  to  those  of  Strich  and  Veillard.  However 
in  our  work,  we  emphasized  the  electronic  structures  of  the  individual  orbitals 
and  the  correlation  of  the  orbitals  of  PF#  with  those  of  the  PF*H,  OPF,,  and 
PFS  molecules. 

The  electronic  structure  of  dinitrogen  tetroxide  and  diboron  tetrafluorlde 
were  determined7  in  SCF  calculations  involving  a moderately-sized  basis  set. 

Since  one  of  these  isoelectronic  molecules  (Na0*)  exhibits  a planar  equilibrium 
structure  in  the  gas  phase  while  the  other  (BaF*)  apparently  has  a rather  small 
rotational  barrier,  their  conformational  stabilities  were  calculated.  Our 
analysis  indicated  that  the  conformational  stabilities  were  determined  primarily 
in  terms  of  lone-pair  interactions  between  nonadjacent  oxygen  or  fluorine  atow 
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along  with  lone-pair  donations  Into  the  central  a * bond.  The  electronic 
structure  of  these  molecules  was  interpreted  in  considerable  detail  and 
it  was  shown  that  previously  suggested*  electronic  configurations  for 
nitrogen  tetroxide  (In  which  the  N0a  units  are  bonded  only  by  it  Interactions) 
are  dissociative  and  hence  of  little  Importance.  Our  conclusions  were  in 
agreement  with  the  available  spectroscopic  data  on  this  molecule. 

Some  of  these  ideas  were  developed  in  a book  entitled  "Electron  Density 
in  Molecules  and  Molecular  Orbitals"  which  was  published  from  our  laboratory 
in  1975,  with  the  proceeds  being  freely  donated  to  Vanderbilt  University. 

In  this  book,  we  have  pointed  out  how  electron-density  plots  may  be  ured 
to  explain  and  interpret  quantum-mechanical  results.  Two  of  the  more  im- 
portant features  of  this  book  are  (a)  the  idea  of  "molecularly  invariant" 
molecular  orbitals,  as  described  in  a preceding  paragraph  and  (b)  the 
effect  of  internal  rotation  on  electron-density  distributions  of  the  various 
valence  orbitals  in  a molecule. 

B.  Calculations  with  the  Core  Electrons  Omitted.  A great  deal  of 
effort  was  spent  in  this  area  and  much  of  it  was  not  published,  since  we 
were  unhappy  with  the  results  and  saw  no  reason  to  clutter  up  the  literature 
with  it.  Many  of  the  original  ideas  (such  as  the  lumped  core)  suggested  in  the 
Grant  Application  simply  did  not  work  out  in  the  long  run,  although  for  a 
while  they  looked  promising  enough  for  us  to  sink  considerable  effort  into 
them.  Much  the  same  thing  could  be  said  about  our  attempts  to  improve  the 
various  semi -empirical  techniques  so  as  to  make  them  less  empirical  and  more 
generally  applicable.  Of  this  kind  of  work,  only  one  investigation  was  published. 
This  paper  presents  a description  of  the  new  DVMO  (directed- valence  molecular 
orbital)  method;  which  is  based  on  the  idea  that,  for  the  results  of  an 
approximate  MO  method  to  be  invariant  upon  rotation  of  the  molecular  coordinate 


system,  there  is  no  requirement  that  directed- valence  effects  (e.g.  the 
difference  between  o and  ir  p-type  orbitals)  be  neglected  in  the  one-center 


Fock-matrix  elements.  This  new  method,  in  which  these  effects  are  accounted 


for,  was  applied  to  a group  of  diatomics  and  found  to  give  reasonably  good 
results.  However,  even  though  this  was  the  best  of  our  endeavors  along 
this  line,  we  do  not  particularly  recommend  our  new  method  for  general 
use  or  further  development  since  it  doesn't  lead  to  significant 


When  Dr.  Patrick  Coffey  joined  our  group  in  May  of  1974,  we  were  feeling 


very  desperate  about  what  route  to  follow  in  order  imitate  or  obviate  the 


core  electrons  of  a SCF  calculation.  We  had  considered  the  pseudopotential 


and  effective-potential  techniques,  but  only  in  passing  because  the 


available  literature  seemed  to  indicate  that  this  approach  did  not  afford 


an  appropriate  basis  for  the  development  of  a generally  applicable  practical 
method  for  carrying  out  LCAO-MO-SCF  calculations.  However  as  a last 


desperate  gamble.  Dr.  Coffey  looked  into  this  method  to  see  if  it  could 


be  developed  in  such  a way  that  it  would  be  suitable  for  routine  quantum' 


chemical  studies.  Fortunately,  he  found  a way  to  do  this,  an  approach  which 


has  been  described  in  the  first  paper  from  our  laboratory  dealing  with  our 


new  NOCOR  (Neglect  of  Core  ORbitals)  method.  In  this  first  paper 


method  was  applied  to  the  main-group  atoms,  ranging  from  atomic  number  2 


through  36,  as  well  as  to  the  molecules  Ca,  Sia,  Gea,  and  PFa.  The 


findings  on  these  atoms  and  molecules  were  compared  to  conventional  SCF 


results  in  comparable  basis  sets  and  to  experiment  for  the  Ca  and  PFa 


molecules.  We  were  particularly  pleased  that  our  electron-density  plots 


of  the  valence  orbitals  of  the  diatomic  molecules  (Ca,  Sla,  and  Gea) 


showed  that  the  spacial  distributions  outside  of  the  core  region  appeared 


to  be  very  similar  to  those  obtained  from  the  equivalent  full-SCF  calculations 
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The  next  paper12  on  the  NOCOR  method  dealt  with  the  application  of  this 
approach  to  the  four  PX*  phosphorus  halides  (X  » F,  Cl,  Br,  and  I),  with 
the  geometries  being  optimized  and  the  force  constants,  orbital  energies, 
Mulliken  populations  and  dipole  moments  being  calculated  using  a minimum- 
Slater  valence  basis  set. 

These  first  two  papers  represented  a practical  test  of  our  ideas  of 
how  pseudopotential  theory  might  be  applied  to  everyday  quantum-chemical 
calculations;  and  these  two  publications  showed  that  the  approach  was 
worthy  of  further  consideration.  A problem  that  often  comes  up  in 
theoretical  work  is  the  question  of  how  much  effort  should  be  put  into 
theoretical  developments  before  the  underlying  ideas  have  been  tested  as 
to  their  usefulness  and  practicability.  In  our  laboratory,  we  usually 
take  the  general  approach  of  first  testing  the  usefulness  of  a theoretical 
advance  while  the  theory  is  still  in  a state  of  flux  and  is  far  from  being 
worked  out  well.  If  the  preliminary  practical  tests  show  merit,  we  then 
attempt  to  investigate  the  theoretical  details. 

Our  first  careful  theoretical  Investigation  of  the  new  computational 
procedure  (the  NOCOR  method)  was  described  In  a paper19  which  appeared  in 
1975.  In  this,  the  theory  of  pseudopotentials  and  effective  potentials  was 
Investigated  and  properly  extended  to  cover  the  approach  we  had  developed. 

This  theoretical  study  indicates  the  following:  (a)  increasing  the  size  of 

the  basis  set  requires  a concomitant  increase  in  the  complexity  of  the  model 
potential  used  to  describe  the  core.  A minimum-Slater  basis  set  works  well 
in  these  calculations  because  the  atomic  and  molecular  eigenfunctions  are 
"locked"  into  essentially  the  same  shape.  However,  the  use  of  a more  flexible 
basis  set  (such  as  a double-zeta  set) , requires  a more  accurately  defined  shape 
of  the  model  potential  in  the  valence  region,  and  in  this  study  we  found  that 


at  least  two  mathematical  terms  (instead  of  a single  one  for  the  single* 
zeta  case)  were  required  for  this  purpose,  (b)  We  did  not  find  any 
evidence  for  failure  of  the  froz«n-core  approximation,  whereby  the 
core  model -potential  of  an  atom  (in  the  form  of  an  atom  or  ion,  or  present 
in  a molecule)  is  considered  to  be  unchanged  by  the  differing  environment  due 
to  the  various  distributions  of  the  valence  electrons  in  the  different 
systems.  This  conclusion  indicates  that  there  is  no  need  for  core -polarization 
potentials,  as  has  been  suggested  by  others.  Furthermore,  (c)  this  work 
has  indicated  that  the  NOCOR  procedure  is  about  as  accurate  as  the  corresponding 
conventional  SCF  calculation  in  determining  molecular  structures,  vibrational 
spectra,  dipole  moments,  Mul liken  populations,  and  the  influence  of  d 
orbitals. 

In  the  next  theoretical  paper, Xh  which  is  still  in  press,  a 
new  expression  for  the  Phillips-Kleinman  pseudopotential  is  proposed  for  many- 
valence-electron  molecules.  According  to  this  approach,  each  valence 
molecular  orbital  is  treated  in  terms  of  its  orthogonality  constraints  with 
respect  to  the  core  orbitals.  This  is  in  contrast  to  our  previous  approach 
and  those  of  all  others,  in  which  an  approximate  average  pseudopotential 
based  on  the  valence  atomic  orbitals  of  the  free  atoms  is  employed.  Our 
new  formulation  of  the  problem  results  in  an  appreciable  increase  in 
precision,  particularly  for  the  valence -orbital  energies.  In  addition, 
we  derived  general  expressions  for  the  local  terms  in  the  effective 
potential,  which  are  not  restricted  to  any  special  analytical  form.  We 
have  now  shown  that,  if  the  core  is  described  precisely  in  terms  of  this 
local  potential  function  representing  only  the  Coulombic  and  exchange  inter- 
actions and  if  the  correct  Phillips-Kleinman  pseudopotential  is  derived 
from  our  new  formulation,  the  entire  effective  potential  is  well  accounted 
for.  In  view  of  this  work,  we  now  believe  that  the  general  approach 
represented  by  the  NOCOR  method  has  been  demonstrated  to  be  fundamentally 
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sound  and  ought  to  be  generally  applicable  for  virtually  all  kinds  of  SCF 
calculations.  In  this  paper,  the  new  approach  has  been  tested  on  the  I* 
molecule;  and  the  results  obtained  with  it  are  compared  with  our  previous 
approaches  and  with  the  results  fron  an  equivalent  full  SCF  calculation. 
The  manuscript  of  this  important  paper  is  presented  as  Appendix  I in 
this  Final  Report.  It  now  seems  quite  clear  on  the  basis  of  this  paper 
that  the  proposals  for  further  quantum- chemical  work  outlined  in  the  Grant 
Application  for  AF0SR-77-314S  (under  which  we  are  now  operating)  are 
well  founded  and  should  be  suitable  for  practical  development. 

In  a series  of  three  papers**"1*  which  appeared  in  1976,  we  have 
applied  the  NOCOR  method  to  practical  problems  of  calculating  molecular 
wavefunctions . In  one  of  these  investigations,**  a calculation  was 
carried  out  on  the  P*  and  Pa  molecules  and  their  Interconversion.  In 
this  work,  the  molecular  geometry  was  optimised  and  the  resulting  bond 
lengths  were  found  to  be  close  to  those  measured  experimentally.  Upon 
allowing  d functions,  it  was  found  that  the  electronic  hybridisation  of 
each  phosphorus  atom  in  the  P*  molecule  changes  from  s 1 ' *°p*' 10  to 
S* * 77p* * °*d° * 19  resulting  in  considerable  polymerisation  of  the  valence- 
electronic  distribution,  with  the  total  electronic  charge  shifting 
away  from  the  atoms  towards  the  center  of  the  P*  tetrahedron  and  also  in 
the  P-P  axial  regions.  Unfortunately,  due  to  an  inappropriate  setting 
of  the  lower  limit  for  calculating  integrals,  a mistake  was  made  in  the 
molecular-orbital  correlation  chart  relating  the  P«  molecule  to  the  pair 
of  Pa  molecules  resulting  from  pulling  a P*  molecule  apart  along  its  Ca 
axis.  This  error  has  been  corrected  in  an  erratum.** 

The  next  paper17  in  this  group  deals  with  NOCOR  calculations  on  the 
tetracarbonyls  of  nickel,  palladium  and  platlmm.  We  felt  that  it  was 
very  important  to  carry  out  soon  a study  on  transition-metal,  compounds. 
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since  gone  persons  in  this  field  think  that  any  pseudopotential  method 
simply  cannot  achieve  reasonable  results  when  applied  to  the  compounds  of 
the  transition  metals.  In  this  work,  we  have  compared  the  full-SCF  results 
on  tetracaxbonylnickel  with  the  NOCOR  calculations  on  this  molecule  and 
found  good  agreement.  However,  we  found  that  the  outermost  orbitals  of 
the  palladium  and  platinum  carbonyls  were  in  reversed  order  with  respect 
to  those  found  for  the  nickel  compound  and  we  spent  considerable  space 
in  the  paper  arguing  how  this  might  have  come  about.  However,  a French 
scientist  has  now  carried  out  a full-SCF  calculation1*  for  tetracarbonyl- 
palladium  and  has  shown  that  the  ordering  of  the  molecular  orbitals  with 
respect  to  energy  in  this  compound  is  the  same  as  in  its  nickel  analog! 

This  scientist  has  therefore  entitled  his  manuscript  "Failure  of  a 
Pseudopotential  Calculation  for  Pd(CO)*"  in  reference  to  our  work.  It 
turns  out  that,  during  the  time  we  were  engaged  in  this  calculation,  a 
minor  eTTor  was  made  in  the  mathematical  description  of  the  core  d orbitals, 
which  are  found  in  the  core  palladium  and  platinum,  but  not  nickel.  There- 
fore the  nickel  calculations  turned  out  correctly  but  there  were  some 
minor  mistakes  which  were  big  enough  to  change  the  ordering  of  the  two  outer- 
most valence  orbitals  in  the  calculations  on  the  palladium  and  platinum 
analogs.  We  are  presently  in  the  process  of  repeating  these  calculations 
and  have  already  found  for  palladium  that,  when  the  core  d orbitals  are 
represented  by  correct  mathematical  expressions,  the  so-called  "failure" 
of  the  pseudopotential  method  disappears.  We  are  happy  that  the  error  was 
a human  one  and  is  not  an  inherent  disadvantage  of  the  NOCOR  method. 

tt 

The  third  paper  of  this  group  deals  with  the  use  of  the  pseudopotential 
method  in  calculating  barriers  to  Internal  rotation.  In  this  work,  the 

i 

barrier  to  internal  rotation  as  computed  by  the  NOCOR  method  for  ethane, 
methylsllane,  methyl germane,  and  methylstannane,  were  compared  with 
conventional  SCF  results  for  all  but  CH,SnH,  and  in  each  case  with  experiment. 

j 
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As  expected,  the  NOCOR  results  turned  out  well.  Analysis  of  the  NOCOR  barrier 
(based  solely  on  the  valence  orbitals)  Indicates  that  the  pseudopotentlal 
formulation  (which  is  Included  in  the  one-electron  part  of  the  Fock  matrix) 
effectively  Incorporates  the  entire  core/valence  interaction,  Including 
the  omitted  two-electron  contributions.  Therefore,  the  pseudopotential 
method  has  now  been  shown  to  work  for  calculating  internal-rotational 
barriers;  and  this  occurs  because  of  the  relative  invariance  of  the  omitted 
core  and  core-core  electronic  contributions  to  the  internal -rotation  process. 
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I THE  STRUCTURE  OF  MOLECULAR  EFFECTIVE  POTENTIALS  IN  COMPOUNDS 

OF  HEAVY  ELEMENTS, WITH  APPLICATION  TO  I2 

Carl  S.  Ewig,  Roman  Osman,  and  John  R.  Van  Uazer 

Department  of  Cheristry,  Vanderbilt  University,  Nashville,  Tennessee  37235 
(Received  ) 

A new  expression  of  the  Phlllips-Klelnman  pseudopotential  for  many— 
valence -electron  molecules  is  proposed.  All  remaining  terms  in  the  valence 
Hamiltonian,  represented  by  a local  potential,  are  evaluated  from  ab  initio 
expressions  based  on  the  form  of  the  atomic  Fock  operator.  Sample  calculations 
are  reported  for  the  orbital  energies,  equilibrium  bond  length  and  vibrational 
frequency  of  the  ground  state  of  the  Ia  molecule. 

I.  Introduction 

The  concept  of  effective  potentials  and  pseudopotentials  for  studying 
atomic  and  molecular  valence-only  electronic  structure  has  in  recent  years 
been  shown  to  be  valid  for  an  extremely  wide  variety  of  different  properties 
and  systems.  In  a previous  paper*  we  showed  that  in  particular  the  one-electron 
form  of  the  Phillips-Kleinman  pseudopotential*' * may  be  rigorously  applied  to  self- 
consistent— field  calculations  of  the  valence  electronic  structure  of  many— 
valence-electron  atoms  and  molecules.  The  primary  advantages  of  this  approach 
are  that  (1)  the  exact  core-valence  orthogonality  of  each  molecular  orbital 
may  be  represented  by  the  pseudopotential,  rather  than  being  Implicitly  main- 
tained at  a value  representing  the  ground-state  neutral  atom,  as  some  other 
valence-electron  approaches  require, *'*»* and  (2)  the  remaining  terns  In  the 
effective  Hamiltonian  are  greatly  simplified,  being  well  described  by  f rations 
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which  are  monotcnically  decreasing  In  the  valence  region  and  very  nearly  the 
same  for  all  symmetries  of  valence  orbitals.  That  this  latter  condition  is 
not  generally  true  of  other  techniques  haa  been  clearly  shown  by  the  work  of 
Kahn  and  Goddard,1  ir-tar  alia. 

The  purpose  of  the  present  paper  is  to  resolve  some  ambiguities  in 
our  previous  formulation  of  the  Phllllps-Kleinman  expression  for  molecules. 
Second,  we  will  suggest  some  more  rigorous  and  accurate  expressions  for  the 
local  "model"  potential  representing  the  core-valence  Coulomb  and  exchange 
operators.  This  potential,  together  with  the  pseudopotential,  we  will  refer  to 
"effective  potential'.'  The  description  of  the  local  potential  is  of  particular 
concern  due  to  the  observation1  that  simple  analytical  expressions  for  the 
model  potential  seem  to  be  Inappropriate  for  use  with  double-seta  and  larger 
valence  basis  sets.  Some  representative  valence-electron  results  for  the 
orbital  energies,  equilibrium  bond  length  and  vibrational  frequency  of  la  in 
a minimal  basis  set  are  reported  and  cocpared-to  the  results  of  an  equivalent 
all-electron  SCF  calculation. 


II.  The  Molecular  Phllllps-Kleinman  Expression 

The  requirement  that  a molecular  orbital  be  normalised  and  orthogonal 

to  soma  set  of  functions,  usually  designated  as  the  core  orbitals,  may  be 

shown1 ***to  be  energetically  equivalent  to  relaxing  the  orthogonality 

constraint,  but  adding  to  the  Hamiltonian  a generalised  pseudopotential  term, 
GPP 

V . Defining  first  the  molecular  projection  operator 


E ' 


i 


I ■ 

B. 

I 

E- 


A nuclei  orbitals  i 1 

p " 1 1 I ♦*><♦<*! 

k i CX  CK 


(1) 


where  1*  the  i**1  core  function  on  nucleus  k,  the  one-electron  form  of 


..GPP 


for  Hartree-Fock  calculations  on  systems  with  one  valence  orbital1  is  then 
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where  F is  the  appropriate  Focfc  operator  and  ey  is  the  eigenvalue  of  F corresponding 
to  the  valence  orbital.  Further,  since  it  is  Assumed  that  the  A*k  are  chosen  to 
approximate  the  true  core  solutions  of  F,  so  that  FA*k  - eck^ck'  Eq’  ^ reduces 
to  the  more  familiar  one-electron  form 


pp  nuclei  orbitals  . 

V = I £ U .>(e 

k i ck  y 


(3) 


If  a valence  wavefunction 


satisfies  the  secular  equation 


0ii 

FA  « c A 

v v v 

than  Eq.  (4)  may  be  written  in  the  equivalent  form 

<F  + VPX  " VV 


<4) 


(5) 


Her*  X^»  the  pseudowavef unc t ion , may  have  a much  simpler  fractional  form  ►h— 

This  follows  from  the  fact  that  the  presence  of  V®*  or  VPP  in  the 

secular  equation  causes  all  the  A^  to  be  degenerate  with  x^,  so  that  xj 

is  arbitrary  in  the  core  region  within  a linear  combination  of  the  A* 

ck 

Hence  it  may  be  chosen  to  be  of  a convenient  form  such  as  a single  Slater 
orbital,  or  otherwise  shaped  in  the  core  region.  Here  we  will  simply  take 
X*  to  be  the  outermost  Slater  function  used  to  describe  the  atomic  A* 
each  nucleus.  This  provides  a computationally  convenient  definition  of 
in  the  core  region,  and  resolves  the  ambiguity  with  respect  to  the 
contributions  of  the  A*k»  With  this  type  of  constraint,  will  not  be 
identical  to  the  exact  valence  wavefunction,  but  the  valence  energy  will 


18 


still  converge  properly.  Thus  computation  of  the  complex  nodal  structure 
of  4*  in  the  core  region,  required  by  conventional  atomic  end  molecular 
SCF  calculations,  may  be  obviated. 

For  free  atoms  and  molecules  with  only  one  valence  orbital  (or  one  of  each 

syamatry),  the  above  procedure  is  straightforward.  However  in  other  molecules 

it  becomes  ambiguous,  since  there  Is  then  a different  e*  for  each  x^«  Other 

authors,  using  expressions  analogous  to  (5)  have  used  an  arbitrary  "naan" 

value7  of  ey,  or  implicitly  kept  it  constant  at  its  value  in  the  free  atom.*""* 

In  our  previous  work,  we  have  set  cy  equal  to  the  lowest  of  the  molecular 

orbital  energies,  found  iteratively. 

PP 

Ideally  V should  be  formulated  such  that  Eq.  (5)  is  satisfied  axactly  for 

each  This  may  be  done  formally,  using  the  fact  that  the  x*  are  normalised 

* 

and  orthogonal,  by  writing  the  molecular  pseudopotential  as 


*\fpp 


valence 
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l*  i 


£ t 
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|>»  (6) 


where  is  the  i**1  valence  eigenvalue  of  F.  designates  this 

more  general  form  of  the  pseudopotential , which  Is  a function  of  the  complete 
set  of  valence-orbital  energies.  In  the  iterative  solution  of  the  self- 
consis tent-field  equations,  the  quantities  ej  and  may  be  conveniently 
chosen  as  the  eigenvalues  and  eigenfunctions,  respectively,  of  the  preceeding 
iteration.  At  convergence,  Eq.  (5)  is  then  satisfied  for  each  valence 
molecular  orbital,  e*. 

It  should  be  noted  that  the  operator  defined  by  Eq.  (6)  is  nonhermit ian,  a 
condition  which  presents  some  difficulties  in  its  direct  implementation . However 
inspection  of  this  form  ofV*P  shows  that  in  the  molecular-orbital  basis  it 


•rnememm* 
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contains  off-dlagonal  matrix  elements  coupling  only  valence  orbitals  of  the 
same  symmetry  type,  and  that  these  elements  are  generally  at  least  two  orders 
of  magnitude  smaller  than  the  diagonal  terms.  Therefore  for  the  present  study 
we  have  constructed as  given  by  Eq.  (6)  In  the  molecular-orbital  basis,  and 
then  neglected  the  off-dlagonal  terms.  Sample  computations  using  this 
expression  are  given  below  in  Section  IV. 


III.  Construction  of  Atomic  Model  Potentials 


The  Fock  operator  In  Eq.  (5)  depends  on  the  coordinates  of  all  the 
electrons  in  the  system,  and  Is  hence  per  ae  inappropriate  for  valence- 
electron  calculations.  A fundamental  postulate  used  in  such  calculations 
is  the  existence  of  at  least  one  function  W1  such  that,  if  Eq.  (5)  is 
•atisfied,  then  for  each  valence  orbital  of  a given  atom 


Cr  + v"  + «*>)£  - ££  <7) 

1 * 
where  xy  i»  the  atomic  psaudowavef unction  and  F now  depends  explicitly  only 

on  the  coordinates  of  the  valence  electrons.  (Note  that,  for  free  atoms, "If 

may  be  chosen  equal  to  V**  for  each  atomic  symmetry  type.)  Since  Eqs.  (4)  and 

(5)  represent  the  conditions  that  the  total  energy  be  an  extremum  with  respect 

to  arbitrary  variations  of  4*  and  xj  respectively,  Eq.  (7)  is  equivalent  to 

requiring  that  the  energy  of  the  valence  electrons  also  be  an  extremum  with 

respect  to  variations  of  the  x^*  Thus  subsequent  molecular  calculations  using 

this  model  potential  are  limited  to  basis  sets  composed  of  functions  which 

precisely  satisfy  the  atomic  secular  equation.  In  a molecular  calculation,  W 

la  taken  to  represent  the  superposition  of  the  potentials  of  the  constituent 


atoms.  In  our  earlier  work,  we  have  used  for  each  atom  a single  expontlal 
function  of  the  form 
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Wexp(r)  - N(l-«"ar)/, 


where  r Is  the  distance  from  the  nucleus  and  N is  the  nunfcer  of  core 
electrons,  a nay  be  determined  by  optioizing  the  values  in  atomic  SCF 
calculations.  We  now  ask  what  is  the  most  general  form  of  W1. 

One  approach  would  be  to  exploit  the  nodeless  properties  of  each  x^  to 
Invert  Eq.  (7),  i.e. 

* ■ H - i - /P)*K-  » 


1 


1 


Since  all  the  quantities  on  the  right  side  of  Eq.  (9)  are  known  from  atomic  SCF 
results,  this  generates  a local  potential  that  may  be  subsequently  employed  in 
molecular  computations.  However  it  is  easy  to  see  that  this  will  be  successful 


only  at  the  Hartree-Fock  limit,  since  using  small  analytical  basis  sets, 

Eqs.  (5)  and  (7)  need  not  be  satisfied  at  all  points  in  space.  In  particular, 
they  need  not  be  satisfied  in  the  outer  or  valence  region  of  the  atom,  which  is 
the  most  important  for  bond  formation.  Trial  calculations  analogous  to  those 
described  below  verified  that,  at  least  for  minimal  basis  sets,  Eq.  (9)  gives 
only  a relatively  crude  molecular  potential  function.  This  and  similar 
problems  that  may  arise  from  using  expressions  analogous  to  (9)  have  also  been 
discussed  by  Melius  and  Goddard." 

A more  general  condition  for  defining  W*  is  to  require  operator  equivalence. 


that  is 


(F  + V^P  + W1)^  - (F  + VPP)X^. 


(10a) 


or,  equivalently. 


<F  + Wi)X,1,  - fA 


(10b) 


f 

J 


21 


Note  that  Gq.  (10)  Is  equivalent  to  requiring  that  an  arbitrary  variation  in 


Xy  cause  the  same  change  in  both  the  total  and  the  valance-electron  energies 
However  this  energy  change  need  not  be  zero,  and  hence  this  need  not  be 
the  solution  of  an  atomic  secular  equation. 

To  calculate  we  first  define  explicitly  the  atonic  Fock  operator  F 
for  closed  shells  as* 


-l5A  - - + 

^ r 


closed 

£ 


£ <2J  -K  ) + f Z <2J  -K  ) 

j 3 3 j 3 3 


where  A is  the  Laplacian,  Z the  nuclear  charge,  and  f is  the  fractional  occupancy 
of  the  open  shell,  if  any.  For  open  shells. 


_ closed  open 

F - -JsA  - - + £ (2J  -K  ) + / £ (2aJ  -bK  ) 

r j 3 3 j 3 3 


where  a and  b are  constants  for  a given  electronic  configuration.  Siailarly, 


F for  closed  shells  is  given  by 

closed 

* „ valence  . _ 


open 

valence 


* -%A  - 7 + £ (2J.-K.)  + / £ (2J.-K.) 

r j J 3 j J 3 


and  for  opan  shells 


closed 

valence  A . 


open 

valence 


F - -JsA  -7  + £ (2J.-K  ) + / £ (2aJ  -bK  ). 

J 3 3 J 33 


We  also  define  J and  K , as  well  as  J and  K , to  be  the  Coulomb  and  exchange 

V V V V*  ° 

operators  defined  over  the  and  xy  respectively,  so  that 


Jv4  - 
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K xj 
vAv 


~ 1 
jv4 


xV 

''if  ”lT 


(16) 

(17) 

(18) 


Further  let  Jc,  J°  etc.  be  the  operators  over  closed  and  open  valence  orbitals 
v v 

respectively . 

* 

Noting  that  the  only  difference  between  F and  F lies  in  the  two-electron 
operators,  it  is  convenient  to  divide  the  W*  functions into  two  terms.  The  first 
is  the  total  interaction  between  the  core  and  valence  electrons; 


■ [«->c  - 


(19) 


where  J and  K are  defined  to  include  all  the  core  orbitals  of  the 
c c 


summations  in  (11)  and  (12).  The  second  term,  W*,  is  due  to  valence  electrons 
in  the  core  region.  Solving  Eq.  (10)  for  W^,  using  Eqs.  (11)-(14)  we  have  for 
closed  valence  orbitals , 

< * K - K - K>-  K - + 4.  ^,^4 


(20) 


and  for  open  valence  orbitals. 


wj  -j[2(J^  + ap;  ~ Kv  - bA°+  K®+  bjT^jxJ/**  (21) 


assuming  all  the  core  orbitals  are  closed.  The  total  potential  Wl  is  then 
W1 


simply  W*  + W*  for  a given  x* 


In  forming  the  molecular  effective  potential,  it  should  be  noted  that 


m - 
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Che  pseudopocentlal  defined  by  Eq.  (6)  is  equally  applicable  Co  both  atonic 
and  molecular  pseudowavefunctions.  In  a molecule  such  as  I, , ic  Is  useful  to 
choose  Che  valence  atcnic  orbitals  of  Che  free  atom  Co  be  Chose  with  Che 
maximum  principal  quantum  number  (i.e.  5 for  Iodine),  and  Che  valence  orbitals 
of  the  nclecule  to  be  those  dominated  by  these  atomic  valence  orbitals.  The 
and  of  the  nolecule  are  given,  to  an  excellent  approximation,  by  the 

eigenvalues  of  the  closed-shell  Fock  matrix  for  the  free  atoms  in  the  double- 
Hamiltonian  formalism  of  Roothaan.*  As  in  our  earlier  work,1*10  take  the  model 
potential  for  molecular  calculations  to  be  the  superposition  of  the  model  potentials 
of  the  constituent  atomic  cores.  Finally,  the  Interatomic  core-core  interactions 
are  represented  by  a reduction  in  each  nuclear  charge  by  N,  the  number  of  core 
electrons,  in  computing  nuclear-nuclear  repulsion  energies. 

IV.  Sample  Computations 

To  illustrate  the  use  of  the  above  expressions  for  both  the  pseudopotential 
and  local  model  potential  in  determining  molecular  properties,  we  have  performed 
both  all-electron  and  effective  potential  minimal-basis-set  calculations  on  the 
ground  state  of  the  I2  molecule.  This  structure  was  chosen  since  it  is  well- 
known  compound  of  reasonably  heavy  elements  and  is  also  amenable  to  full  SCF 
calculations. 

All  results  reported  here  were  obtained  using  an  atom-optimized  set  of  Slater 

orbitals.  Each  core  Slater  function  was  expanded  into  three  gaussian  functions 

\ 

and  each  valence  orbital  was  expanded  into  four  gaussians.  The  a parameter 
in  Eq.  (8)  was  obtained  by  fitting  the  atomic  ground-state  orbital  energies  as 
has  been  described  previously,1*10*11  and  for  iodine  has  the  value  3.1814  a.u.-1 
The  properties  chosen  for  comparison  are  the  orbital  energies  of  the  five  occupied 
valence  orbitals,  the  equilibrium  bond  length  and  the  vibrational  frequency. 

The  effective-potential  calculations  in  each  case  required  less  computer  time 
than  the  corresponding  full  SCF  calculation  by  a factor  of  150. 
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The  second  and  third  columns  of  Table  ! show  the  effect  of  using  a pseudo- 
potential  of  the  form  given  by  Eq.  (3),  setting  ty  to  be  the  lowest  molecular 
valance-orbital  energy  e compared  to  Eq.  (6)  in  which  all  the  e*  are 

found  iteratively.  The  most  apparent  difference  is  seen  to  be  an  improvement 
in  the  accuracy  of  the  orbital  energies  relative  to  the  full  SCF  results. 

In  computing  the  local  potential  functions,  it  is  clear  from  Eqs.  (19).  (20) 

and  (21)  chat  there  will  be  a different  W*  for  each  symmetry  of  valence  basis 

function.  For  example,  in  the  minimal  basis  of  the  iodine  atom  there  will  be 

one  for  x*  and  a different  one  for  These  hopefully  should  be  smoothly 

v v 

decreasing  in  the  valence  region,  and  W8  should  be  nearly  equal  to  W^.  This  is 
indeed  seen  to  be  the  case  from  Table  II,  which  shows,  in  the  first  four  columns 
of  data,  the  core  and  valence  potentials  and  Wy  appropriate  to  the  x^  and  x^ 
atomic  pseudowavef unctions. 

Although  W8  and  could  be  used  separately  in  a molecular  valence- 
electron  calculation,  we  will  use  here  a single  function,  W,  to  represent  both 
in  order  to  draw  a valid  comparison  with  the  use  of  the  exponential  approximation, 
wexp.  Thus  we  will  simply  define  W as  the  least-squares  fit  to  W*(x^) * for  all 
atomic  symmetries,  i,  at  each  point  in  space,  so  that  for  the  iodine  atom 


w8^)4  + WP(tf) 
<xj>4  + <x£)4 


This  function  is  shown  in  the  fifth  column  of  Table  II,  and  Wex**  is  presented 
for  comparison  in  the  sixth  column. 

In  performing  molecular  computations  in  a contracted  gaussian  basis  set, 
it  is  clearly  desirable  to  have  the  local  potential  expressed  as  an  expansion 
in  gaussian  functions. 


W(r)  * N[l  - Ic1exp(-e1r*)]/r 


(23) 
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This  expansion  was  found  to  converge,  after  inclusion  of  five  or  more  terms, 
to  approximately  three  significant  figures  of  the  properties  given  in  Table  I. 
The  weighting  factor  employed  was  the  denomimator  in  Eq.  (22)  multiplied  by  r*. 
The  coefficients  and  exponents  of  this  expansion  are  shown  in  Table  III. 

The  properties  found  using  this  potential  with  the  pseudopotential  defined 
by  Eq.  (6)  (using  the  molecular  pseudowavef unctions)  are  shown  in  the  fourth 
column  of  Table  I.  The  differences  from  the  results  obtained  with  W*3^  are 
again  seen  to  be  relatively  small. 

V.  Discussion 

The  advantages  of  employing  the  Phllllps-Klelnman  expression  in  molecular 
calculations  were  stated  in  Section  I.  Further,  its  use  leads  to  Eqs.  (5)  and 
(10).  These  equations  in  turn  allow  construction  of  the  single  and  general 
functions  of  the  form  (19),  (20)  and  (21).  These  local  potential  functions 
have  the  properties  that  (1)  the  functional  form  of  W1  is  not  constrained  to 
to  any  analytical  approximation,  and  (2)  the  atomic  and  functions  need 
not  be  solutions  of  an  atomic  secular  equation,  thereby  permitlng  the  use 

of,  for  example,  molecular ly-optimized  basis  sets. 

The  purpose  of  the  present  paper  has  been  to  describe  the  origin  of 
some  of  the  terms  appearing  in  the  molecular  effective  potential.  Although 
the  results  reported  here  are  somewhat  more  precise  than  we  had  previously 
thought  possible  for  this  row  of  the  Periodic  Table,  it  seems  very  likely 
that  more  precise  valence-electron  calculations  are  possible.  Parameterization 
schemes  *» * * have  been  found  to  work  well  for  compounds  of  light  elements, 
but  become  extremely  complex  for  the  heavier  elements  since  they  must 
implicitly  simulate  Eq.  (3).  However  the  accuracy  of  the  valence-electron 
approximation  itself  is  not  yet  well  established  for  compounds  such  as  Ia. 

Other  pseudopotential  results  for  Ia  are  not  yet  readily  available,  although 
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Kahn  et  al.**  report  an  equilibrium  bond  length  of  2.81  A using  a somewhat 

o 

core  complex  basis  set  in  an  MCSCF  formalism.  The  experi rectal  value  is  2.667  A.1 
Recently  Barthelat,  Durand  and  Serafini,**  using  a doubie-zata  basis  set, 
have  obtained  orbital  energies  for  I*  differing  from  the  corresponding  full 
SCF  results  by  an  average  of  12.51,  and  an  equilibrium  bond  length  and 
vibrational  frequency  differing  by  82  and  17.42  respectively. 

An  apparent  minor  lmprovmant  that  might  be  made  in  the  computational 
approach  described  here  would  be  the  use  of  a different  V*  function  for  each 
atomic  symmetry  type.  The  problems  that  arise  In  this  type  of  calculation,  due 
to  the  fact  that  a basis  function  on  nucleus  A need  not  have  a simple  atomic 
symmetry  in  the  region  of  nucleus  B,  have  been  treated  approximately  by 


several  authors.* **,ia  However,  since  the  W®  function  is  very  similar  to 
the  \P  function.  It  seems  likely  that  their  use  in  molecule;  calculations 
would  be  quite  insensitive  to  the  exact  fashion  in  which  this  angular- 
momentum  dependence  were  treated. 
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